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We introduce the method of the double pole QCD sum rules and we study the p(2S), tf>(2S) and 
T(2S) mesons. We get the masses and decay constants of the mesons and show a prediction for the 
decay constant of the p(2S) meson given by (177 ± 13) MeV. 
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I. INTRODUCTION 

In 1979, Shifman, Vainshtein and Zakharov (SVZ) [lj created the successful method of QCD sum rules 
(QCDSR), which is widely used nowadays and they won the Sakurai prize in 1999 for this discovery. With 
this method, we can calculate many hadron parameters like: hadron's mass, decay constant, coupling 
constant and form factors in terms of the QCD parameters like: quark masses, the strong coupling and 
non-perturbative parameters like quark condensate and gluon condensate. The main point of this method 
is that the quantum numbers and quarks hadron's content are represented by an interpolating current, 
where we get a correlation function as a time order operation of this current. The determination of the 
ground state of hadron's mass, we use the two point correlation function, where this correlation function 
has two point of view, the QCD side and the phenomenological side. On the QCD side, the correlation 
function can be written in terms of a dispersion relation, where the spectral density doesn't have poles, 
on the other hand, the phenomenological side has a sum of the all states that coupling with the current. 
To represent the phenomenological side as a form of a dispersion relation the SVZ method uses an Ansatz 
that the spectral density can be represented by the form "pole + continuum" . 

For the p meson's spectrum the proposal of SVZ is a good approach, due the large decay width of the 
p{2S) or p(1450) and the possible existence of a light p(3S) or p(1570) with decay width about 144 MeV, 
Fig.(jl|, that allow to approximate the excited states of this theory as a continuum, on the other hand, 
the J/ip meson's spectrum is very different and a proposed type "pole+ pole+ continuum" seems more 
realistic, Fig. dJ). 

Alternative proposals for the phenomenological spectral density have been used in the QCDSR to study 
the case where the hadron ground state has a large decay width [2] , where the pole is replaced by Breit- 
Wigner function. The case where the excited states are considered in QCDSR there are two methods: 
the Maximum Entropy Method Q and Gaussian Sum Rule with "pole+ pole+ continuum" Ansatz Q. 
In Maximum Entropy Method there are studies for the p meson nucleon J/ip meson Q and if 
meson Q. In Gaussian Sum Rule there is a study for mixed states of the glueballs and scalar mesons. 

In lattice QCD the study of excited states is a recent area pL where there are a lot of studies devoted 
to 7r(25) meson [§], p meson excited states @|, charmonium pol-[l2|. nucleon excited states p^ - [l5| and 
exotic charmonium spectrum [lfj. In addition, the excited states have been studied recently by several 
approaches like: QCD's Bethe-Salpeter Equation [llj for tt(2S) and p(2S), light-front quark model fl8l . [l9l ] 
for p(2S), r] c (2S), ip(2S) and the bottomonium analogous. The ip(2S) has been studied in QCDSR as a 
hybrid meson [2(| using the "pole -I- continuum" Ansatz. 

There are a lot of motivation to study the excited states that belong the charmonium spectrum. States 
recently discovered like Y(4260) and Y(4660) are a great ex amp le of this lack, considering theories that 
Y(4260) has been proposed as a bound state of J/ip — fO [2l| and Y(4660) has been interpreted as a 
bound state of J/tp(2S) - /0, [HSH show that Y(4660) is an excited state of Y(4260). Another point 
is that Z + (4430) would be an excited state of X + (3872) and Z^(10610) would be an excited state of 
X+ (10100) [13. 
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The simple method "pole+ polc+ continuum" Ansatz was used in lattice QCD [15| and the authors 
have shown a problem in which the nucleon current coupling more with excited state than its ground 
state. In Gaussian Sum Rule the "pole+ polc+ continuum" Ansatz [1] was applied and they get the 
masses of the mixed states of the glueballs and scalar mesons. 

In this paper, we study the excited state using the "pole+ pole+ continuum" Ansatz in QCD sum 
rules and we apply in three cases: the p(2S), ip(2S) and T(2S) mesons and we calculate their masses and 
decay constants. 
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FIG. 1: (Left) The radial excited states of the p meson [2a [26|]. The p(1540) and p(1900) are omitted from 
PDG's summary table, but p(1900) is a good candidate for p(3S) was predicted by Refs. [2a. [27j| . For /p(1540), 
its existence is not predicted by usual qq model. (Right) The radial excited states of the J/i/j meson [25l. \WL [281] . 



II. THE METHOD 

To implement our method, we first use the invariant part of the generic correlation function and we 
apply the Borel transformation: 



oo 

B[U phen (q 2 )} = J dsp phen (s) 



(1) 



where r = 1/M 2 . We consider the following spectral density to the phenomenological side: 



p ph ™(s) = XlS(s - ml) + X 2 2 S(s - m 2 2 ) + p c ° nt (s)6(s - s' Q ), 



(2) 



where m± is the mass for the ground state and m-i is the mass for the first excited state and s' a mark 
the begin of the continuum states. 

Inserting Eq.© in Eq. ([!]), we get the expression: 



n phen {r) = \\e- m * T + \\e-< T + / dsp phen (s) 



(3) 



where W hen (r) = B[W hen {q 2 )]. 

Using the quark hadron duality, we get the double pole QCD sum Rule, 



U QCD {1 



where, 



(4) 



U QCD ( T ) = J ds e~ ST p OPE (s) + higher order condensates, 



(5) 
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where s™ m is a QCD parameter. 

As usually done in QCDSR to obtain the mass of the hadron, we take the derivative of Eq. (|4]) with 
respect to r and we get the new equation: 

- m\\\e- m ^ - m\\\e- m ^ = -^U^ cd (t). (6) 

(XT 

We can observe that the equations Eq. (gj) and Eq. ^ they form a equation's system in the variables, 

A(t) = \le~ m *\ (7) 



B{t) = \ 2 2 e- m ^ T . (8) 

Solving the equation's system Eq. (@} and Eq. © writings in terms of the functions A{t) and B(t), 
we easily get: 

A(T) = —2 2 ' I 9 ) 



jn^(r)+nQ^(r)mj 

B{T) = 2 2 ' l 1U ) 

mj — 77^2 

where: DF{t) = £F(t). 

To eliminate the dependence of the Ai coupling that appears in Eq. ([5]), we take a derivative of this 
equation with respect of r and divide the result by Eq. ([9]). The result of this procedure is given by the 
Eq. (fTTj) . To eliminate A2 coupling the procedure is analogous and the result is given by the Eq. (fl"2|) . 



/ DHQ CD (t) m| + DDW CD (r) 
mi -V DUQ CD (T)+UQ CD {T)m 2 2 ' ( ' 



j DDUQ CD (r) + DUQ CD (t) mf 
TO2_ Y nO cc (r)mf + DnQ CD (r) ' ( ^ 

In the first view the Eq. lfTT]) and Eq.([T2l) suggest a system of non-linear equations for the masses mi 
and m2, that could be extracted in independent way. On the other hand, using Eq. (|lll) to obtain an 
mi expression, it reproduces the same result given in Eq. (|12[) . It is a very interesting result, that is 
impossible to isolate the masses mi and m2- In that way, our sum rule the mass of the ground state is a 
function of the mass of the excited state, where we use in our analysis the equation Eq .([TT|). 

III. SUBTRACTION OF THE CORRELATION FUNCTION 

For the qq vector mesons, on the QCD side the correlation function has a form: 

U^q) =i fdf-x e*»(0| T{j M (x)jt(0)} |0> = (q»q„ - q'g^U^q 2 ) , (13) 



where j^{x) = q a {x) 7^ q a (x) . 

On the phenomenological side we use: 



(0|i/.(0)|V(g)> = A.mve^(g), (14) 
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where /y is the decay constant and my is the hadron's mass. Inserting Eq. (fT4")) in Eq. fTB")) , we get: 



P 

H/j.v(q) = (l/jiQu — m v9^v) — — — — 2 + excited states contribution. (15) 



Comparing the Eq. flTSl) with Eq. (fT5")) we notice that they have different structures. We choose the g^ v 
structure that it doesn't depend on the choice of a gauge on the phenomenological side. We can write 
the invariant part of the correlatores of the Eqs. (fT3")) and (fT5")) . in a form: 



n(, 2 )=/ds^, (16) 

S — q l 



Considering the subtracted correlator [2S 



% 2 ) = i% 2 ) - n(o), (17) 

we get: 

oo 

flfo 8 )^ / 'd8-fi°Lr, (18) 

J s(s - q z ) 



s and 

For the sum rules of p, J/tp and T mesons we use the Rcf.|2S 



where we use for the sum rule the function, M| ' . Comparing the p(s) with spectral density given by 
Eqs.© and ©, we get: p OPE (s) = p OPE (s)/s and p phen {s) = 5 phen (s)/s. 



IV. RESULTS 



In this work we use the following parameters: a s = 0.3, m q — 4MeV, m c = 1.3 GeV, mb = 4.5 GeV, 
(qq) = — (0.23) 3 GeV 3 and (g 2 G 2 ) = 0.88 GeV 4 . For the study of the meson's mass, is used a curve in 
two dimensions where the stability as a Borel mass is studied for three level curves: Mi, solid line, that 
it supplies a double pole that contributes with 90% of Eq.©, M2, dash line, that contributes with 60% 
and M3, dash-point line, that contributes with 40%. 

This imposition that the pole contribution defined the Borel window is a recent progress of the QCDSR. 
The use of this criterion for the pentaquarks 9 + (1540) and S (1836), has been shown that the QCDSR 
not support a pentaquark state 30]. Recently the same criterion was used in Ref.(22| and they have 
shown that is impossible to interpret Y(4260) as a hadron molecule J/ip — Jq. 



A. p(2S) and p Sum Rule 

Using the data of p meson's spectrum Fig. dU) , the s' Q = (mz + A') 2 parameter has A' = 0.4 GeV. But 
this value supply unstable values for the decay constant. The best value for A' is 0.1 GeV. In Fig. we 
show in the level curves of the mass Eq. (jlip . where p meson's mass, mi as a function of the p(2S) mass, 
m,2. We notice that with the increase of the Borel mass the curve of the mass converge to the curve with 
Borel mass M3. Seemingly, none point seems privileged, but the mass difference, m,2 — mi, it is only 
possible experimental value A = 0.69 GeV for = 1.46 GeV. In this case, we have 7711(7712) = 0.77 GeV. 
So for this sum rule it is necessary to do some estimate for the separation of the states 2S and IS, to do 
the prediction of the mass. 

For the calculation of the decay constant, we use the experimental values mi = 0.77 GeV, 7772 = 1.46 
GeV and A' = Q.lGeV. In Fig. ([3]), we show the decay constant of the p meson, solid line and the decay 
constant of the p(2S) meson, dash line. We see the value for the p meson decay constant has a plateau 
on value 200 MeV. This value is in agreement with the experimental value, given by 216 MeV [25|]. For 
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FIG. 2: The p mass as a function of the p(2S) mass for A' = 0.1 GeV. 



the p(2S) meson decay constant has a plateau with a value of 165 MeV. Considering uncertainty with 
respect to A' parameter due large width of the p(2S) and uncertainty of the p(3S) mass, we varying this 
parameter in the range: A' = (0.15 ± 0.05) GeV, we have a prediction for its decay constant given by, 



/ p(2S) = (177±13)MeV. 
For the p meson decay constant, we get: f p — (200 ± 2) MeV. 



(19) 




FIG. 3: The decay constant of the p meson, solid line, and p(2S), dash line, as a function of the Borel mass with 
A' = 0.1 GeV. 



In Fig.Q we show the pole contribution as a function of the Borel mass, comparing the sum rule with 
two poles for s' — (1.46 + 0.1) 2 GeV 2 , dash line, with the sum rule of usual one pole sum rule for the 
so = (0.77 + 0.7) 2 GeV 2 , solid line. We see that the sum rule with two poles has a Borel window larger 
than the usual case. 

In Fig.© we show the relative contribution of the p meson, dash line, and p(2S), solid line, that 
compose the double pole. We see that for low values for the Borel mass the contribution of the meson p 
it reaches 90%, but for Borel mass M 3 it reaches a lowest contribution of the 67%. 



B. ip(2S) and t/j Sum Rule 

In the study of the J/ip, we use the data of J/ifj meson's spectrum Fig. (J}, where A' = 0.3 GeV. In 
this case, we get the mass of the J/ip for M3 Borel mass, given by mi(3.7) = 2.9 GeV, Fig. ©. 

On the other hand, the calculation of the decay constant, we use the experimental values m\ =3.1 GeV, 
rri2 = 3.7 GeV and A' = 0.3 GeV in Eqs.© and ([TU|) . In Fig. (J7J) we see the values for the decay constant 
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FIG. 4: The pole contribution as a function 
of the Borel mass, comparing the two poles 

, dash line, and one pole sum rule, solid FIG - 5: The relative contribution of the p 
j me meson, dash line, and p(2S), solid line. 
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FIG. 6: The J/ip mass as a function of the ip(2S) mass. 



of the J j ip meson, solid line, has a good stability for the Borel mass above 3 GeV and its value is close 
the its experimental value, that is about 416 MeV. For the decay constant of the ^f(2S) meson, dash line, 
has the same stability window and its value is (320 ± 10) MeV. Considering A' parameter in the range: 
A' = (0.3 ± 0.05) GeV, we get: 

U(2S) = (327 ±37) MeV. (20) 

This result is in agreement with the experimental value, given by (297 ± 3) MeV. 

In Fig.© we show the contribution of the pole as a function of the Borel mass, comparing the sum rule 
with two poles for s' — (3.7 + 0.3) 2 GeV 2 , dash line with the usual sum rule for s = (3.1 + 0.56) 2 GeV 2 , 
solid line. We see that the sum rule with two poles has a Borel window more larger than the usual case. 
In Fig.© we show the relative contribution of the J/4> meson, trace line, and ip(2S), solid line, that 
compose the double pole. We see that for low values for the mass of Borel the contribution of the meson 
J/ijj it reaches 95% for the Borel mass M\. 

C. T(2S) and T Sum Rule 

In the study of the T, we use A' = 0.33 GeV. In this case, we get the mass of the T(IS') for A/3 Borel 
mass, given by mi (10.02) = 9.24 GeV, Fig. (fT0|). 

On the other hand, the calculation of the decay constant, we use the experimental values mi = 
9.46 GeV, m 2 = 10.02 GeV and A' = 0.33 GeV in Eqs.© and (TTUl). In Fig.(HH) we see the values for the 
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FIG. 7: The decay constant of the J/ij) meson, solid line, and ip{2S), dash line, as a function of the Borel mass 
with A' = 0.3 GeV. 
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FIG. 8: The pole contribution as a func- 
tion of the Borel mass, comparing the two FIG. 9: The relative contribution of the 
poles, dash line, and one pole sum rule, J/t/j meson, dash line, and tf>(2S), solid 
solid line. line. 
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FIG. 10: The T mass as a function of the T(2S) mass. 



decay constant of the T meson, solid line, has a good stability for the Borel mass above 5 GeV and its 
value is close the its experimental value, that is 714 MeV. For the decay constant of the T(2S) meson, 
dash line, has the same stability window and its value is, 



h(2S) = (530 ± 30) MeV. 



(21) 
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This result is in agreement with the experimental value, given by (497 ±3) MeV. In this case the 
uncertainty with respect to A' parameter does not exist, because the T(3S) is a narrow ressonance, with 
full width of 20 KeV. 
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FIG. 11: The decay constant of the T meson, solid line, and T(2S), dash line, as a function of the Borel mass 
with A' = 0.33 GeV. 



In Fig.(fT2j) we show the contribution of the pole as a function of the Borel mass, comparing the sum rule 
with two poles for s' = (10.02+0.33) 2 GeV 2 , dash line with the usual sum rule for s = (9.46+0.56) 2 GeV 2 , 
solid line. We see that the sum rule with two poles has a Borel window more larger than the usual case. In 
Fig. flTBl we show the relative contribution of the T meson, dash line, and T(2S), solid line, that compose 
the double pole. We see that for low values for the mass of Borel the contribution of the meson T(2S) it 
reaches 92% for the Borel mass M\. 




FIG. 12: The pole contribution as a func- 
tion of the Borel mass, comparing the two 

poles, dash line, and one pole sum rule, FIG - 13: The relative contribution of the 
solid line. f meson, dash line, and T(2S), solid line. 
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V. CONCLUSIONS 



In this work we have presented a new method to QCD sum rule with double pole and we have applied 
this method for three problems well established in QCD sum rules, the study of the p, J/tp and T mesons. 

We have shown that the obtaining of the hadron's masses supply close values to the experimental 
masses and we have shown a prediction for the decay constant of the p(2S) and we have obtained the 
decay constants of the ip(2S) and T(2S) and their ground states that are shown in Table fl] 



TABLE I: Decay Constants of the 2S states and IS states in MeV. 
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Ref.[10j 


Ref.[25j 
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200 ± 2 216.37 


268 






216 ±5 


p(2S) 


177 ±13 128 


155 








J/1> 


416 ±11 






399 ±4 


416 ± 5 


i>{2S) 


327 ± 37 




371 


143 ± 81 


295 ±3 


T 


695 ± 11 




546.6 




714 ±5 


T(25) 


530 ± 30 




583.2 




497 ±3 
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